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Abstract. We calculate numerical solutions and analytic approximations for the intermediate-type 
spectral distortions. Detection of a //-type distortion (saturated comptonization) in the CMB will 
constrain the time of energy injection to be at a redshift 2 x 10^ > z > 2 x 10^, while a detection 
of a j-type distortion (minimal comptonization) will mean that there was heating of CMB at redshift 
z < 1.5x10^*. We point out that the partially comptonized spectral distortions, generated in the redshift 
range 1.5 x 10^ < z < 2 x 10^, are much richer in information than the pure y and //-type distortions. 
The spectrum created during this period is intermediate between y and //-type distortions and depends 
sensitively on the redshift of energy injection. These intermediate-type distortions cannot be mim- 
icked by a mixture of y and //-type distortions at all frequencies and vice versa. The measurement of 
these intermediate-type CMB spectral distortions has the possibility to constrain precisely not only 
the amount of energy release in the early Universe but also the mechanism, for example, particle 
annihilation and Silk damping can be distinguished from particle decay. The intermediate-type dis- 
tortion templates and software code using these templates to calculate the CMB spectral distortions 
for user-defined energy injection rate are made publicly available. 
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1 Introduction 

Cosmic background explorer (COBE/FIRAS) [1] measurements show that cosmic microwave back- 
ground (CMB) follows Planck spectrum to a high precision between 1 < ;c < 1 1, where x = hv/k^T is 
the dimensionless frequency, h is Planck's constant, is Boltzmann's constant and T is the temper- 
ature of the CMB blackbody spectrum. The precision is quantified by the 2cr limits on the chemical 
potential [2] jU < 9 x 10"^ and y-type distortion [3] y < 1.5 x 10"^. Technologically an improvement 
of more than two orders of magnitude over COBE/FIRAS has been possible for some time [4] and 
proposed future experiment Primordial Inflation Explorer (PIXIE) [5] will be able to measure ay-type 
distortion of y = 10~^ or = 5 x 10"^ at 5cr, a more than three orders of magnitude improvement 
over COBE/FIRAS. 

In case there is energy or photon production at a redshift z > 2 x 10^, the photon production and 
destruction through bremsstrahlung and double Compton scattering along with the redistribution of 
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Figure 1. Important epochs related to the creation and evolution of the CMB spectrum. 



photons in energy via Compton scattering on thermal electrons can establish full thermal equilibrium 
[2, 6, 7] and we get a blackbody spectrum with a higher temperature. The redshift z « 2x 10^ therefore 
defines the boundary of the blackbody photosphere. This happens, for example, when electrons and 
positrons, having higher initial entropy and energy density than photons, annihilate at z ~ 10^ ^ 10^. 
Also, when there is energy release during primordial nucleosynthesis, the photon spectrum quickly 
thermalizes. For redshifts between 2xl0^<z<2xl0^, any heating of CMB gives rise to a Bose- 
Einstein spectrum or a //-type distortion, where // is the chemical potential. For redshifts z < 1 .5 x lO'^, 
the Compton redistribution of photons over frequencies is too weak to establish the equilibrium Bose- 
Einstein spectrum and we get a y-type spectrum. For jU-type (3'-type) distortions, if detected, we can 
only put a lower (upper) limit to the time of energy release. However, if there is heating of CMB in 
the redshift range 1.5 x lO'* < z < 2 x 10^, the spectrum depends sensitively on the time of energy 
injection and it is thus possible to put a much more precise constraint on the time of energy release. 
Important epochs in the history of the early Universe, from point of view of the CMB spectrum, ai^e 
depicted schematically in Fig. 1. 

We first define the spectral distortion as a pure redistribution of photons of a reference black- 
body, which is slightly different for /u type distortions from the conventional definition. This definition 
makes it possible to uniquely define the zero point/crossing frequency of the spectral distortion and 
use it to determine the redshift of energy injection. We review the .y-type solution and also calculate 
the regime of validity of this solution in the Appendices. We study the comptonization of y-type 
distortion and calculate the evolution of zero point/crossing frequency of the distortion. As a first 
example, we discuss the shape of the spectral distortions arising from the dissipation of sound waves 
in the early Universe [8-14] and how they can be used to constrain the shape of the primordial power 
spectrum on small scales with wavenumbers 8 < /c < 378 Mpc"'. Additional significant sources of 
heating/cooling in the early Universe include dark matter (WIMP) annihilation and adiabatic cooling 
of baryons [12, see also 7, 15]. Thej-type distortions are also created throughout the late-time history 
of the universe, in particular from the heating of CMB by hotter electrons in the intergalactic medium 
during and after reionization, and by hot intracluster electrons in the clusters of galaxies. The early 
Universe contributions to the j-type distortions, before recombination, are indistinguishable from the 
late-time contributions; y-type distortions are therefore not very useful in constraining possible new 
physics, which can heat the CMB before reionization. As long as the distortions are small, the 3'-type 
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distortions from different epochs and the fi and the intermediate-type distortions just add linearly to 
give the total distortion to the CMB. Evolution of spectral distortions in the early Universe was first 
considered by Zeldovich and Sunyaev [2, 3] and analytic and numerical solutions were computed 
by [16]. These and later calculations [6, 12, 17-19], which included double Compton scattering and 
considered low baryon density Universes such as ours, were focused on /u and y-type distortions, 
although some authors computed the full evolution of the spectrum, including intermediate-type con- 
tributions. We should clarify that in this paper we are always in the non-relativistic regime and are 
not concerned about the relativistic corrections to either the Kompaneets equation [20] or the j-type 
distortions, which have been studied in some detail by various authors [21-32] in applications to the 
very hot (Jg > 10 KeV) gas inside clusters of galaxies. 

In the following calculations we will use WMAP cosmological parameters of ACDM cosmol- 
ogy [33] with CMB temperature today Tcmb - 2.725K, number of neutrinos Ny = 3.046 [34], helium 
mass fraction yne = 0.24, Hubble constant Hq = 70.2km/s/Mpc, matter density n„j = 0.275, baryon 
density Q.h = 0.0458 and zero curvature. 

2 Possible non-standard sources of energy injection in the early Universe 

Detection of energy release in the early Universe is an important source of information about new 
physics. There are several possible theoretical sources from high energy theories. For example, 
in super-symmetric theories and Kaluza-Klein theories it is possible that dark matter was initially 
produced as a long lived next to Ughtest particle in the dark sector and then decays later to the lightest 
particle which acts as dark matter today [35]. An example is a neutralino (B) decaying to gravitino 
(G) in super-symmetric theories. Neutralino in this case is a WIMP which decays into gravitino 
which has only gravitational interaction is thus super weakly interacting or SWIM? [36]. Sterile 
neutrinos can also be a significant component of dark matter and their decay can be a source of 
energy and photons [35]. Other sources of energy release include evaporating primordial black holes 
[37], decaying cosmic strings and other topological defects, cosmic string wakes [38] and oscillating 
super-conducting cosmic strings [39, 40], and small-scale primordial magnetic fields [41]. 

In most of the above examples the decay of an initial particle results in high energy electromag- 
netic and hadronic showers. Initial standard model particles, that ai^e produced as a result of these 
showers, interact with electrons, ions and photons through a rich variety of energetic processes like 
Compton and inverse Compton scattering, pair production, photon-photon interaction etc. (see for 
example [42]) depositing most of the initial energy in the form of heat in plasma very quickly. Some 
of the energy is lost to neutrinos and energetic particles in the windows of low optical depth. This 
energy deposition in turn gives rise to a j-type distortion which then comptonizes. As a second ex- 
ample, we consider decay of an unstable particle as a source of energy injection and show how the 
intermediate-type spectral distortions can constrain the Ufe time of the particle as well as distinguish 
it from energy injection due to a process for which the rate of energy injection as a function of redshift 
is not exponential but a power law, such as annihilation of particles or dissipation of sound waves. 

3 Definition of CMB reference temperature and spectral distortion 

For any given isotropic unpolarized photon spectrum , «(v), where n is the occupation number as a 
function of frequency y, we can define a unique reference temperature. At redshifts z < 2 x 10^, the 
energy exchange between the plasma and radiation is very fast, but the production (and absorption) of 
photons is extremely slow. As a result, the number density of photons does not change due to energy 
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release and the ratio of photon to baryon number density ?iy /?ib is constant to high precision. ' We can 
calculate the total number density of photons in the spectrum 

r 87ry2 

A^= I —^n{v)dv, (3.1) 

where c is the speed of light. We can now define a reference temperature, T, as the temperature 
of the blackbody spectrum «pi which has the same number density N. Thus T = (A'^/^r)'^^, where 

= ^^'^c^//^^^ ' ^ Riemann zeta function with ^(3) « 1.20206. We can now also define 

dimensionless frequency, x = hvl{k^T) and write the total spectrum as 

n{x) - npi{x) + An{x), (3.2) 

where An is the distortion from the reference blackbody with the property that it represents a redis- 
tribution of photons in the reference blackbody, «pi = l/(e^ - 1). The total number of photons in the 
distortion vanishes, j dx x^An = 0. This formalism is especially useful for discussing the distortions 
created in the early Universe solely through the action of comptonization, which just redistributes the 
photons already existing in a previous spectrum, and thus does not change the reference temperature 
defined above. Of course, many other definitions of spectral distortions aie possible by defining a 
different reference temperature [12, 13]. The above definition is just an extension of the conventional 
definition of a v-type distortion to distortion of any shape. 

In the rest of the paper we will use the above definitions of reference temperature and spectral 
distortions. We note that for the case of a Bose-Einstein spectrum «be> the above definition gives a 

different reference temperature than the one given by the usual definition ?ibe = 1/ ^e^BE """^ - 1 j. Let 
us remind that hbe is the equilibrium solution of the Kompaneets equation [20] when Tg = Tbe, where 
Tg is the electron temperature. Our definition corresponds to taking an initial blackbody spectrum 
with temperature equal to the reference temperature T and add to it small amount of energy (keeping 
the photon number constant), which then fully comptonizes creating a Bose-Einstein spectrum with 
a chemical potential ju fully defined by the amount of energy release. It is easy to calculate the 
final temperature of the resulting spectrum using relations for number density and energy density 
of a Bose-Einstein spectrum in the limit of small chemical potential [44] and it is given by ? = 
(Tbe - T)/T = 0.456/1. Thus the spectrum written in terms of dimensionless frequency x = hvlk^T 
is given by, 

1 

eVm — 1 

1 

gA-0.456/i.r+/.i _ \ 

giving the zero point, where the coefficient of n in above equation vanishes making the spectrum 
identical to the blackbody spectrum at temperature T, at xq - 2.19, v - 124 GHz, compared to 
xq = 3.83, V = 217 GHz for a j-type distortion [3], 



xe' 

ny{x) 



(e^ - 1)2 



(3.4) 



'There is an exception, of course, when the energy release mechanism also adds photons, for example Silk damping 
[13,43]. 
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Finally we note that part of the energy release from sound wave dissipation in the early Universe 
goes into blackbody part of the spectrum [13]. The blackbody part of the energy release is easily ab- 
sorbed in the reference temperature T and can be ignored completely in our method, where we define 
the spectral distortions as pure redistribution of photons of a reference blackbody.^ We have plot- 
ted in Fig2. 3'-type and yu-type (Bose-Einstein) distortions having same number and energy density. 
The fractional difference in the effective temperature is plotted, which is defined by the following 
equations. 



n(x) 



1 



1 



Ar _ Teff - r ^ 1 - e"-^ An 
T T X ripi 



(3.5) 



Fig. 3 shows the difference in intensity from blackbody, Aly = ly'^ - ly^ oz x' An, where /^'-^ is 
the intensity of /i or j distorted spectrum, and /y' is the intensity of Planck spectrum at reference 
temperature T defined above. The two type of distortions shown (ju and y) have the same energy 
density. 



4 Creation of y-type distortion, role of Compton yy parameter and validity of the 
solution 

A source of energy injection, for example dark matter decay, will in general lead to a shower of 
particles which will quickly deposit most of their energy in the plasma, and result in an increase in 
the electron temperature as long as the source is on. We will first review the j-type distortions created 
as a result of a source of energy that turns on for a very short time and also discussed the regime of 
validity of this solution. 

4.1 Compton parameter jy and the Compton distortion parameter y. 

The equation describing the evolution of photon spectrum through Compton scattering is the Kom- 
paneets equation [20]. We will work with the dimensionless frequency x and the photon occupation 
number n{x) which is given in the case of blackbody spectrum at temperature T = 2.725(1 -i- z) K by 
fipiix) = I lie" - 1). In the case of a blackbody or a Bose-Einstein spectrum at a different tempera- 
ture, there will be factors associated with the rescaling of temperature. Working with n{x), which is 
invariant with respect to the adiabatic expansion of photon gas, allows us to avoid the extra terms in 
the equations associated with the expansion of the Universe. To simplify equations further, we will 
use as our time variable the parameter jy (not to be confused with the ^'-type distortion parameter y) 
defined as follows: 



yy{Z,Zm^yd = - dz 777; (4.1) 

A™x m,c H{\+z) 



^If we add of order e fractional energy to a blackbody of initial temperature T with part of it going into blackbody 
and part into fi distortion, then for the new temperature T' we have, {T' - T)/T ~ ^ ~ e. Ignoring the energy addition 
to the blackbody part we will use the frequency variable x, while the frequency variable with respect to the new reference 
temperature is x' = hv/(kgT'). The error introduced in the calculation of distortions is thus second order in e, An(x)/n(x) = 
An(x')/n(x') + 0(£'), and can be ignored for small distortions. In addition, this change in blackbody temperature takes no 
part in comptonization, since it also changes the electron temperature by the same amount (Eq. (5.1)). 
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Figure 2. y-type and //-type distortions created by addition of energy £ s bJLjE = 4x10 ^toa blackbody 
with reference temperature T . Fractional difference in effective temperature defined in Eq. 3.5 is plotted. 



where z^ax is the redshift where we start the evolution of the spectrum or the energy injection redshift. 
It is convenient to use radiation temperature T and not electi'on temperature in the definition of x 
and yy, as it makes x independent of the expansion of the Universe after electron-positron annihilation, 
with both r, y oc (1 -Fz). The electron temperature on the other hand evolves in a non-trivial way after 
baryons thermally decouple from radiation at z ~ 500. The total yy parameter for an energy injection 
redshift of Zinj is then given by j'yCzinj) = 3'y(0,Zmj)- The Compton pai'ameter jyCzinj) is plotted in 
Fig. 4 and it can be seen that the contribution to the integral from z < 1.5 x 10^ becomes very 
small, with yy < 0.01. We will present our results as functions of yy and they can be converted into 
functions of redshifts using Eq. 4. 1 and Fig. 4. During radiation domination, we can calculate the 

yy parameter analytically and is given by 3'y(z,Zmax) « ^'m!/ ^""^2^^In'/^''^'' "*" ^^^"''^ ~ "*" ^ 
4.88 X 10"^' [^(1 -I- Zmax)^ - (1 + z)^] > where nw and ?i4He are the number densities of hydrogen and 
helium nuclei today respectively, and Qr is the radiation energy density today in units of critical 
density. Similarly, it is easy to find analytic formulae during the matter dominated era, z ^ 3200, but 
before recombination and also after reionization, when the electron density is again simply equal to 
sum of hydrogen and helium number densities (assuming singly ionized helium). 
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Figure 3. Same as Fig. 2 but difference in intensity from a blackbody with reference temperature T is plotted. 



We also define the Compton j-type distortion amplitude, y, 



y 



= f\Tdyy = - f 

Jo Jzn, 

J~ kBCTT 



ax 



nieC H{\ + z) 



m+z) 

dz = ye, 



dz 



(4.2) 



where, Aj^ = T^/T - 1, and is the electron temperature. The last line gives the familiar result 
in terms of electron pressure relevant for hot electrons, for example, in the clusters of galaxies. The 
Compton parameter yy is a measure of the degree of comptonization including the effect of all three 
relevant processes: recoil, induced recoil and Doppler broadening (the latter determined by yyT^/T). 
The distortion amplitude y on the other hand is a measure of the amount of heating (or cooling). As 
an example, we show in Fig. 5 the expected post-recombination rate of thermal y-distortion injection 
into the CMB. Before reionization, we get negative 3'-type distortions due to comptonization with 
colder electrons which cool much faster with the expansion of the Universe [12] and can be identified 
as Bose-Einstein condensation [15]. Reionization, in addition to increasing ionization, also heats up 
the intergalactic medium. We have assumed = 10^/(1 + z)^'^ K at z < 3 taking into account 
the contributions from the warm hot intergalactic medium (WHIM) [51, 52] and 10"^ K at z > 3. 
Total y-type distortion from reionization is y^Ar^ ~ {kBTe/meC^)Tri ~ 10"^ x 0.1 = 10"^ while from 
WHIM at smaller redshifts, it is expected to be ~ 10"^ [53-56]. These j-type distortions should be 
detectable by PIXIE [5] while the proposed experiment Cosmic Origins Explorer (COrE) [57] with 
similar resolution but significantly higher sensitivity than PLANCK spacecraft [58] and ground based 
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dz [kBO-jneT] / [mecH{\ + z)]. The 
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drop in the plot at z ~ 1000 is due to the depletion of electrons because of recombination. The total energy 
released is divided approximately equally between y-type, intermediate-type and //-type distortions for the 
case of Silk damping with spectral index of initial power spectrum close to unity [e.g. see 15] and for other 
mechanisms with similar rate of energy injection, e.g. dark matter annihilation. For energy release from 
exponential decay of particles, the division of energy is less democratic and depends sensitively on the lifetime 
of the particle. Recombination [45-47] was calculated using the effective multilevel approach [48] following 
pubUcly available codes HyRec [49] and CosmoRec [50]. 



experiments ACTPol [59] and SPTPol [60] should also be able to detect ^-fluctuations in the WHIM. 
We have also assumed that reionization starts at z ~ 20 and ends at z ~ 8. Uncertainties in the 
details of reionization and temperature evolution of the inteigalactic medium make it impossible, at 
least at present, to disentangle the pre-recombination and post-recombination contributions to the y- 
type distortions. There will also be contributions from the second order Doppler effect from baryon 
peculiar velocities [61, 62], not shown in the figure, before during and after reionization. These 
contributions during and after recombination are calculated in [13]. 

With the above definitions, we can now use the standard form of Kompaneets equation [20] (tak- 
ing care to distinguish between the actual electron/effective photon temperature and the blackbody 
temperature T used to define the variable x): 

-x^\n + n^ + ——]. (4.3) 



dyy dx \ T dx 

The first term in the parenthesis describes the downward scattering of photons due to electron recoil, 
second term is for induced scattering while the last term describes diff"usion of photons in energy due 
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Figures. Post-recombination (rough) estimate of the thermal y-distortion injection rate -(l+z)dy/dzis plotted, 
it is approximately equivalent to the energy release in the redshift interval 5z ~ Z- Before reionization starts, 
the distortions are dominated by cooling of CMB or Bose-Einstein condensation (BEC) due to comptonization 
with colder electrons giving a negative y distortion. During and after reionization the intergalactic medium is 
heated to temperature Te > 10"* giving a much larger positive y distortion of amplitude ~ 10"^. There will also 
be similar magnitude positive contributions from the second order Doppler effect (not shown above) arising 
due to peculiar velocities of electrons. 

to Doppler effect and thus depends on the electron temperature. At high frequencies « «; 1 and the 
induced scattering term can be neglected. If a source of energy raises the election temperature such 
that » r for a very short time, «; 1 , an analytical solution, for the initial condition n{x, yy = 
0) = npiiT) (blackbody spectrum at temperature T), of the Kompaneets equation can be obtained by 
approximating the recoil terms n + with the initial blackbody spectrum npi + Wpi^. The fact that 
the Planck (and in general Bose-Einstein) spectrum is an equilibrium solution of the Kompaneets 
equation ([20], [63], [44]) gives n + ri^ k -dn/dx. The resulting equation can be transformed into heat 
diffusion equation , the analytic solution of which is the well known j-type distortion [3], Eq. (3.4). 
We derive these results in Appendix B in a way which clearly illustrates the regime of validity of the 
solution. 

5 Evolution of j-type distortion 

We will now consider the problem of comptonization of an initial j-type distortion created by an 
energy source which turned on for a very short time, i.e., instantaneous energy injection. This problem 
will illustrate the main physics we want to investigate. We will explore the more realistic cases of 



-9- 



continuous energy injection in the next sections. However, it is possible that such short lived sources 
may actually exist. In fact, one such example can be found in standard cosmology. The decay of 
primordial ^ Be to ''Li lasts for a very short time at z ~ 30000 [64] and gives rise to exactly the type 
of distortions we calculate in this section (although energy released in Be decay is too small to be 
of observational interest). The spectrum we get from ^ Be is the yy = 0.04 spectrum and is given 
approximately by Eq. (5.3) withy ~ 10"'^. 

The y-type spectrum should evolve towards the Bose-Einstein equilibrium solution with time. 
The equilibrium electron temperature or effective photon temperature Te is given by [65, 66] 



The effective temperature for the linear j-type spectrum for energy injection AE/E is given by [44] 
Te = T(l + 5.4y) = T{1 + 1.35AE/E). The exact temperature will be slightly higher if y is not small. 
For a Bose-Einstein spectrum the effective temperature is « T{1 + 0.456//) w T{1 + 0.64AE/E). 
This temperature is established very fast compared to any other relevant process with a characteristic 
time of~ lsatz = 10^, ~ 10 orders of magnitude faster than the expansion rate of the Universe at that 
time. The electrons will thus always be maintained at the effective photon temperature given by Eq. 
5.1. As the y-type spectrum evolves towards the Bose-Einstein spectrum, the electron temperature 
should decrease. 

5.1 Numerical solution in the intermediate era, 1.5 x 10"* < z < 2 x 10^ 

To follow the evolution of the spectral distortions starting with the j-type distortion, we must solve 
Eq. 5.1 and Kompaneets equation Eq. 4.3 simultaneously. Numerically we proceed as follows. We 
take small steps in time y^ using Kompaneets equation with constant Tg given by Eq. 5.1 for the 
spectrum at the beginning of the step. We then calculate the final electron temperature using Eq. 5.1 
and iterate, with linearly decreasing between the initial and final values. We found that a step size 
of 6yy = 0.001 at yy < I and 6yy = 0.01 at > 1 was sufficiently accurate. With our iterative 
procedure the error in energy conservation is < 1% at < 10. Fig. 6 show the cooling of the photon 
spectrum as the Bose-Einstein distribution is approached. Initially the temperature drops rapidly from 
the y-type value of AT^/T = 5Ay and is close to the linear Bose-Einstein value of 2.56y at = 1. 
Fig. 6 shows that for small distortions, An/n «; 1 and \{Ts - T)/T\ «; 1, the equilibrium Bose- 
Einstein spectrum is reached at » 1 and the spectrum is very close to the equilibrium aiyy ~ 1. 
This conclusion does not depend on the amplitude of the distortion, y, at all because the processes 
responsible for comptonization, Doppler broadening and recoil, are defined by the pai^ameters yyT^/T 
and yy respectively. 

Difference in observed intensity from a blackbody is shown in Fig. 7 and Fig. 8 shows the 
fractional difference in the effective temperature with respect to T as the photon distribution moves 
from the .y-type spectrum towards the Bose-Einstein spectrum. An interesting feature is that the 
zero point, defined as xq such that n(xo) = np\(xo), moves from the j-type distortion value of xq = 
3.83 to Bose-Einstein value of xq = 2.19. The maxima and minima of the intensity distortion also 
move towards smaller frequencies as comptonization progresses. Zero point xq is plotted in Fig. 
9. Frequencies of maxima and minima, x^in, ^max are also plotted in Fig. 9 and the corresponding 
intensities, A/min, A/^ax in Fig. 10. The Bose-Einstein spectrum at x > 10 is in fact established very 
quickly. By yy = 0.2 the spectrum is very close to the Bose-Einstein spectrum conxsponding to the 
electron temperature T^iyy) at x > 10. At Vy > 0.2, the spectrum at high x remains Bose-Einstein 
and tracks the electron temperature as the effective radiation/electron temperature decreases. Fitting 
formulae for xq, Xmin, x^ax ai"e given in Appendix A. 
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Figure 6. Evolution of electron temperature for initial spectrum with y = 10"^. Initial fractional difference 
in temperature ATg/T = (T^ - T)/T = 5.4 x 10"^. The final temperature for the Bose-Einstein spectrum is 
also marked at AT^/T = 2.56 x 10"^. By yy - 0.01 the temperature is significantly different from the initial 
temperature and by = 2 it is very close to the Bose-Einstein spectrum temperature. This and subsequent 
plots in this section can be interpreted as snapshots in the evolution of initial y-spectrum starting from initial 
energy injection redshift Zm-dx to redshift z, so that the x-axes is y^ = yyiZmax) - yyiz)- Alternatively, it can be 
interpreted as the final spectrum today resulting from the energy injection at redshift Zinj, so that yy = yyiZmj)- 

5.2 Analytic solution in the weak comptonization limit 

We can find an analytic solution for the evolution of an initial y-type distortion by expanding n(x,yy) 
around y-y = 0, using the Taylor series expansion, Eq. (B.l). The initial spectrum is the ^'-type 
spectrum with amplitude y, n(x, 0) = npi{x) + yny(x), and the initial electron temperature is given 
by the equilibrium temperature, Eq. (5.1), Aj-^ w 5Ay. Substituting the initial spectrum in the 
Kompaneets equation gives us the first term in the Taylor series (assuming y «c 1 and only keeping 
terms linear in y), 



where ny{x) is the j-type spectrum defined in Eq. (3.4), and fyix) is the same function which appeared 
in Eq. (B.3) and is given in the appendix, Eq. (C.l). We thus have the solution for the weakly 
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Figure 7. Intermediate-type spectra: Difference in intensity from a blackbody at temperature T is plotted. 
Bose-Einstein and y-type spectra for same value of energy injection are also shown. Initial spectrum is a pure 
y-type distortion (Eq. (3.4) with y - 10"^) and is labeled jy - 0. The curves in order of increasing (non 
zero) yy correspond to energy injection redshift of Zinj - 1 .56 x lO"*, 3.33 X 10"^, 4.67 X lO'*, 6.55 x lO'*, 1 .03 x 
10^, 1.45 X 10^,2.04 X 10^. The 'zero' point, the frequency x where the intensity equals that of blackbody at 
temperature T moves from the y-type distortion value of 3.83 to Bose-Einstein value of 2.19. 



comptonized spectrum, 



n{x,yy) K ripiix) + y [ny{x) + jy {5Any{x) + /v(jc))] 
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where we have the following simphfied expression and large x and small x limits. 
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(5.3) 



(5.4) 



The correction fy{x), of course, similarly to ny{x), conserves photon number, dxfy{x)x - 0. 

The next terms in the Taylor series can be calculated iteratively by taking successive -derivatives 
of Kompaneets equation. We give the recursion relations to calculate any higher order term in 
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Figures. Intermediate-type spectra: Fractional difference ^ ^ ^~r~«7 in effective temperature (Eq. (3.5)) 
is plotted. Bose-Einstein spectrum for same value of energy injection is also shown. Initial spectrum is a pure 
y-type distortion (Eq. (3.4) with y - 10"^) and is labeled = 0. The curves in order of increasing (non zero) 
correspond to energy injection redshift of z,„; = 1 .56 x 10^ 3.33 x lO'', 4.67 x 10^, 6.55 x 10'', 1 .03 x 10^ 1 .45 x 
10^, 2.04 X 10^. The 'zero' point, maxima and minima of the frequency move towards smaller frequencies as 
comptonization progresses. 



Appendix C. The -derivatives of the electron temperature are also required and are easily cal- 
culated using Eq. (5.1). The first two derivatives are given by, dATjdyy\y^=o w -21.45j, and 
d^Ar^/djy^ly =0 « 323. 6j. Intermediate distortions for case of continuous energy release, for ex- 
ample particle decay/annihilation. Silk damping, are easily obtained from these analytic formulae by 
linearly adding (integrating) the spectra for different yy with appropriate weights. Analytic solution 
including first three terms are quite precise (~ 1% error) for yy < 0.05 deteriorating to ~ 10% errors 
at = 0.1. Numerical and analytic solutions are compared in detail in Appendix C. 



6 Application: Amplitude, slope and shape of the primordial power spectrum on 
small scales 

The solutions given in the previous section would be directly applicable if the energy injection occurs 
over a very short period of time. It is more likely, in reality, that the energy release happens over 
an extended period of time, for example, decay of particles or dissipation of sound waves. The 
final spectrum for continuous energy injection would be a superposition of spectra for all values 
of yy parameter, with appropriate weights decided by the dependence of rate of energy injection on 
redshift, and we must calculate the spectrum for each model of energy release numerically. The shape 
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Figure 9. Evolution of the zero point xo defined by «(jco) - npi(xo), and frequency of minima and maxima of 
the intensity of distortion, Xmin and JCmax- xq can be used to pinpoint the redshift of energy injection in case a 
distortion in CMB spectrum is detected. 



of the power spectrum and the value of xq will still contain information about the energy release as 
a function of time, in addition to the total amount of energy released. This is in contrast to the pure 
//-type (or y-type) distortion which only contains information about the total energy injected. 

One of the most important sources of heating in standard cosmology is the dissipation of sound 
waves in the early Universe because of the shear viscosity (and at late times also due to thermal 
conduction) on small scales. The dissipation of sound waves resulting from thermal conduction was 
first calculated by Silk [8], Peebles and Yu [71] included shear viscosity and Kaiser [72] included 
the effect of photon polarization. The spectral distortions arising from the dissipation of sound waves 
have also been studied by many authors in the past using approximate estimates [9-11] and a precise 
calculation was done recently by [13]. Primordial perturbations excite standing sound waves in the 
early Universe on scales smaller than the sound horizon [73-75]. Diffusion of photons from different 
phases of the waves, of wavelength of the order of diffusion length, gives rise to a local quadrupole 
which is isotropized by Thomson scattering (shear viscosity). The effect on the photon spectrum is 
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Figure 10. Minimum and Maximum of the intensity of distortion from the reference blackbody, A/y 



just the averaging of the blackbodies of different temperature and gives rise to y-type distortions [76]. 
The y-type distortions can then comptonize, fully or partially, giving rise to a yu-type distortion or an 
intermediate-type spectrum [9, 1 1]. We should also mention that the adiabatic cooling of baryons due 
to the expansion of the Universe gives spectral distortions of an amplitude opposite to those given by 
heating [12]. We include this cooling of baryons (or equivalently small difference of electron temper- 
ature from the effective photon temperature given by Eq. (5.1)) in our calculations; this is, however, 
a very small coiTcction to the amount of heating considered below. Also, the low frequency spectrum 
is affected by bremsstrahlung emission/absorption after recombination [12]. In the frequency range 
of interest to us, v > 30 G¥lz,x > 0.5 and for distortions of interest, with amplitude AE/E > 10"^, 
low redshift bremsstrahlung (and double Compton scattering) can be neglected. 

The precise total spectrum resulting from sound wave dissipation was calculated recently by 
[13] including contributions from the yu-type era, intermediate era and the y-type era. Here we con- 
sider the possibility that the pure //-type distortions created at yy > yy^^^ = 2(z > 2 x 10^) and pure 
y-type distortions created at yy < y^^j^ = 0.01(z < 1.6 x 10^) can be subtracted with high precision 
(see Fig. 7). The exact upper/lower limits (and the resulting intermediate spectrum) will of course 
depend on the ability of the experiment to distinguish a pure //-type (y-type) from a yy^nax^rmm^ 
intermediate-type spectrum. The heating rate due to an initial power spectrum with constant scalar 
index on small-scales is given by [13, 43] 
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Figure 11. The power spectra, Eq. (6.3), with different indices «ss on small scales are shown, for constant 
amplitude. Also shown for reference is the best-fit WMAP power spectrum on large scales with «s = 0.96 
[33, 67, 68] with Ly-a forest extending the constraints to smaller scales [69, 70]. The small scale limits on 
power from COBE/FIRAS measurements of CMB spectral distortions are also shown. There is considerable 
freedom in varying the amplitude and the spectral index of the power spectrum within COBE/FIRAS limits. 
We show this by the curve with 10 times the amplitude and extreme value for the spectral index Wss = 2. PIXIE 
[5] is expected to improve COBE/FIRAS constraints by a factor of ~ 2.5 x lO"*. 



where £ = E/py, E is the total energy in photons and py = aRTcuB^i^ + z)'* is the reference photon 
energy density, is the radiation constant, Tcmb - 2.725 K is the CMB temperature today, is the 
damping wavenumber given by [72, 77] 

k^^ ~ J, ^6H{1 + R)n,(Tj \ 1 + /? ^ 15 

where R = 3ph/4py, ph is the baryon energy density. We have defined the power spectrum of initial 
curvature perturbation in comoving gauge ^ as 

An important point to note here is that for «ss - 1> the energy released between redshifts z\ and Z2 is 
proportional to ln[(l + zi)/(l + Z2)] [15] and it is easily seen that the total energy released at z > 1000 
is divided approximately equally between y, /u and intermediate-type distortions. For smaller spectral 
index Wgs < 1> intermediate-type distortions get bigger and bigger share of the total energy released 



(6.2) 
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Figure 12. Distortion spectrum created by dissipation of sound waves in the redshift range 1.56 x 10 < 
z < 2.04 X 10^ corresponding to the jy parameter 0.01 < jy < 2. Difference in intensity from the reference 
blackbody spectrum with the same number density of photons is plotted. Different plots are for different values 
of small scale spectral index Wss of the initial power spectrum normalized so that they all have the same power 
at the pivot point k - A2 Mpc^' as the WMAP best fit power spectrum with «s = 0.96, which is the value of 
diffusion wavenumber at - 0.1, z = 4.67 x 10"*. These are the distortions that would be left over after 
pure jy and ji type distortions are subtracted and probe shape as well as the amplitude of the small scale power 
spectrum in the range % < {k ~ ku) < 378 Mpc"'. 



with 3'-type distortions comparatively enhanced and //-type distortions comparatively suppressed. The 
exact opposite, of course, happens for larger spectral indices > 1- 

Since we want to compare the shape of spectral distortions for different power spectra with 
similar total energy input, we choose the pivot point - ^0(^7 - 0.1, z - 4.67 x IC*) = 42 Mpc"^ 
corresponding to the approximate (geometric) middle of our redshift range for intermediate-type 
spectral distortions. We choose the amplitude - 1.61 x 10~^ to match the power at k - 42 Mpc"^ 
with the WMAP best fit power spectrum with «s - 0.96 [33]. The spectral index, /iss> on small 
scales, ^ > 8 Mpc"' can be very different from the spectral index on large scales, n^, measured 
by WMAP, for example, in the case of a running spectrum. The small-scale power spectra, with 
different indices Wss, are shown in Fig. 11. CMB constrains the primordial power spectrum only on 
large scales [33, 67, 68], k < 0.2 Mpc"^ This range can be extended to ~ 1 Mpc"' using Ly-a forest 
[69] and is consistent with the WMAP best fit power spectrum pai^ameters [70]. The small scales are 
best constrained by COBE/FIRAS data through Umits on y and //-type distortions [1]. We have shown 
these constraints, assuming n^s = 1 and using the fitting formulae in [13] to calculate the y amd // type 
distortions. These constraints ai^e more than 3 orders of magnitude higher than a simple extrapolation 
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of the WMAP best fit power spectrum to small scales. Thus, there is considerable freedom, from an 
observational viewpoint, for the power spectrum on the small scales to be quite different from the 
extrapolation of the WMAP power spectrum to these scales. The constraints get considerable weaker 
behind the blackbody surface at small scales, k > k^iz = 2 x 10^) = 1.15 x 10"^ Mpc"\ because 
of the suppression of the ju distortion by blackbody visibility function [2]. The visibility function is 
dominated by the double Compton process in a low baryon density Universe such as ours [6] and 
is given by w exp [-(z/2 X 10^)^''^j.^. We sketch this weakening of the constraints by multiplying 

the COBE/FIRAS constraint by the inverse of the visibility function exp [(^/1.15 x 10"^ Mpc~^)^^^j, 
using ki) oc (l + z)^^^- The yu-type distortion, of course, only provide integrated constraints on the 
total energy injected in the //-distortion and its separation into individual contributions from different 
epochs in not possible in practice. 

To calculate the spectrum arising from a continuous source of heating, such as dissipation of 
sound waves, we should solve the Kompaneets equation Eq. (4.3) with a source term on the right 
hand side given by 



dn 

dyy 



ld£_d£ 

4 dz dyy 



lly, (6.4) 



where «y is the j-type distortion given by Eq. (3.4) which is created initially and the factor of 1/4 
comes from the relation between the energy injected and the amplitude y of the y-type distortion. 
We show the intermediate-type spectrum resulting from the dissipation of sound waves for several 
different values of the small scale spectral index 0.5 < ?iss < 1-5 in Figs. 12 (intensity). The amplitude 
of the distortion at low and high frequencies is similar since the total energy released is similar in all 
cases. The shape of the spectrum for different values of spectral index is also similar in the Rayleigh- 
Jeans and Wien tails, as expected for the intermediate type spectra. Fig. 7. But the spectra are very 
different and easily distinguishable near the zero crossing, which is at xq = 3.04, v = 173 GHz 
for ?iss = 1- For > 1 there is more power on smaller scales which dissipate earlier moving the 
xq towai^ds lower values (or towards //-type value of xq = 2.19, v = 124 GHz) and for ?iss = 1.5, 
xq = 2.1 A, V =156 GHz. On the other side, for ?iss < 1> the spectrum moves towards the y-type 
value of Xq - 3.83, v = 217 GHz, and for «ss = 0.5 the zero crossing is at xq - 3.32, y = 189 GHz. 
There is of course more information in the full spectrum than just the zero crossing and a sensitive 
experiment should be able to use the full spectrum to tightly constrain more complicated shapes of the 
small scale primordial spectrum than the simple two parameter (amplitude and spectral index) model 
considered here. The zero crossing, and frequencies of minimum and maximum flux difference with 
respect to the reference blackbody can be fitted by the following simple formula as a function of «ss 
for 0.5 < «ss ^ 1.5 at better than 1% accuracy, 

■XO/min/max(«ss) - ^0 + ^l^ss, (6-5) 

where qq - 3.61, ai = -0.588 for xq, = 2.11, ai - -0.466 for Xmin and ao - 6.24, ai = -0.63 for 



7 Application: annihilation and decay of particles 

There are many sources of heating possible in the early Universe from particle physics beyond the 
standard model, as discussed in the introduction. Different sources of energy injection may have 

^Analytic solution including comptonization and both bremsstrahlung and double Compton processes with percent level 
accuracy is given in [7]. 
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Figure 13. Distortion spectrum created by annihilation of WIMP dark matter (multiplied by 400 to bring 
its amplitude to ~ 10"^) and decay of unstable particles with different decay times (zx^tx) - (3.33 x 
lO'*, 661 years), (4.67 x 10^ 339 years), (7.00x 10^ 152 years), (1.03 x 10^ 70.5 years), (1.45 x 10^ 35.7 years); 
the curves are labeled by zx- The shape of the WIMP annihilation spectral distortion is close to the «ss - 1 
spectrum in Fig. 12. Difference in intensity from the reference blackbody with same number density of photons 
is plotted. 



different dependence on redshift and will give rise to different shapes of intermediate-type spectrum. 
We thus have a way of distinguishing between different types of energy injection mechanisms. We 
will illustrate this by considering annihilation of weakly interacting massive particles (WIMP dark 
matter), which has a power law dependence on redshift/time, and decay of unstable particles having 
an exponential dependence on time. 

For the thermally produced WIMP dark matter consisting of self-annihilating Majorana parti- 
cles,^ the energy release due to annihilation is given by, 

d£ _ mmC^nl^ < crv > 
d7 " Py(l+z)H 

6.9xl0-'»(l+z)-i / 10GeV \ 

^J^+{l+Zec,)/i^+Z)\ mm I' 

where we have assumed velocity averaged cross section < crv >« 3 x 10~^^/(ndm^o) cm^^^"^ [78], 
Q.dm is the dark matter density as a function of critical density today, H is the Hubble parameter, 
ho = Hq/IOO = 0.702, fy is the fraction of energy going into heating the plasma, widm is the mass of 



''For Dirac particles the energy release is smaller by a factor of 2. 



dark matter particle, ?idm is the dark matter number density and Zeq ~ 3234 is the redshift of matter 
radiation equality. At z » Zeq we have d£/dz oc 1 + z. For the dissipation of sound waves in Eq. (6.1) 
under same approximation we have d£/dz oc (1 +2)(3nss-5)/2 p^j. _ ^ ^j^j.^, j^^tter annihilation and 
sound wave dissipation have the same redshift dependence and we expect the shape of intermediate- 
type spectrum for these two cases. More precisely, for dark matter annihilation, approximately 30% 
of the energy released at z > 500 appears as yu-type distortion, 37% as intermediate-type distortion 
and the rest of the energy goes to the j^-type distortion. 

For decay of a particle of mass mx, initial comoving number density nxo, and life time tx with 
/y fraction of energy going into heating of the plasma, we have 

d£ _ _ nxomxc^e''^'''' 
" ''aRTcMB'^Hil+z)^Tx 

To get total distortion of ~ 10"^, we choose fyUxomxc^ - 10"^«RrcMB^(l + z.x), where zx 
is the decay redshift corresponding to the lifetime tx, and during radiation domination we have 
Tx ~ 1/(2(1 -I- zx)^HoD.r^^^). The division of released energy into y, fi and intermediate-type distor- 
tions varies quite dramatically with the lifetime of the particle and illustrates nicely how using the 
intermediate-distortions along with ju and 3'-type distortions can help remove degeneracies associated 
with different energy release mechanisms. For zx = 1.45 x 10^, 21% if released energy appeai^s as 
/i-type distortion and 79% as intermediate type distortions. For - 4.67 x 10^,7 x 10'*, almost all 
of the energy, ~ 99%, 97% respectively, goes to intermediate-type distortions with the most of the 
rest in yu-type distortions. For zx = 1-5 x 10^ the division is 42 : 58 between intermediate and ^'-type 
distortions respectively and yu-type distortions get a negligible share. 

The intermediate-type spectral distortion for WIMP annihilation and decay of an unstable par- 
ticles with different lifetimes zx is shown in Fig. 13 (intensity). The WIMP spectral distortion is 
multiplied by 400 to bring it to the same level as the decaying particle signal. WIMP annihilation 
spectrum is similar to n^^ = 1 spectrum in Fig. 12 as expected from their similar- redshift dependence. 
Also, the WIMP annihilation (power law dependence of energy injection on redshift) and the spectra 
for decaying particles (exponential dependence on redshift) with different lifetimes are distinguish- 
able from each other. There is a small degeneracy for zx * 7 x lO'*; the exponential decay in this 
case has same zero crossing as annihilation, and the shapes of two curves ai^e very close. The yu-type 
distortions in the two cases are very different, as discussed above, and break this degeneracy. Never- 
theless, the intermediate-type spectrum has the possibility of measuring the life-time of the decaying 
particle in addition to the total energy injected into the CMB. The frequencies xq, x^i^, -^max can be 
fitted by the formula for dark matter decay at with Z4 = Zx/10^ for 10^ < Zx ^ 2 x 10^ with better 
than 1% precision, 

'"-Omin/man ~ 

ao + a\ ln(z4) + a2 ln^(Z4) + ^3 ln^(Z4) + ^4 ln'^(Z4), (7.3) 

where for xq the fit coefficients are ao - 3.61, a\ - -0.0644,^2 = -0.017, (33 - -0.1365, (34 - 
0.0343, for x^in they ai'e qq = 2.16,ai = -0.07, c?2 = 0.0146,a3 = -0.1165, ^4 = 0.0276, and 
for Xmin we have £?() = 6.25, ai = 0.0158,^2 = -0.157, 03 = -0.1,^4 = 0.0335. For dark matter 
annihilation, the frequencies are xq — 2.96, Xmin = 1-58, x^ax — 

5.52. 

8 Non-degeneracy among Intermediate-type distortions and a mixture of j and yu-type 
distortions 

Since the intermediate type distortions he in-between y and yU type distortions, an important question 
arises: can an intermediate-type contribution to the CMB spectral distortion be mistaken for a com- 
bination of y and ju type distortions and vice versa? The answer is no, these three types of distortions 
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Figure 14. Intermediate-type distortions from particle decay are shown for two different decay times, as in 
Fig. 13. Also shown is the least squares fit of the mixtures of y and fi type distortions (dotted curves) which 
approximate these intermediate-type distortions at PIXIE frequencies. 



are non-degenerate with each other and can be distinguished. This is clear by looking at Fig. 7. 
The intermediate distortions (AT/T) for 0.01 < < 2 are almost constant at both low and high 
frequencies. The j-type distortions, on the other hand, rise at high frequencies as AT/T oc x. The 
magnitude of the jU-type distortions similarly increases at low frequencies with AT/T oc 1/x. Any 
mixture of pure y and yu type distortions will thus have much greater slopes than the intermediate type 
distortions, and in principle they can be separated from each other. 

We show in Fig. 14 the same spectrum as in Fig. 13, for the energy injection due to particle 
decay for two different decay times, 2x = 1.45x10^,3.33x10^. Also shown is a combination of ju and 
y-type distortions, which approximate these intermediate-type spectral distortions in the least squares 
sense. Thus, a sensitive experiment should be able to distinguish between /i, y and intermediate-type 
distortions when the distortions are detected at high significance. In particular, it should be possible 
to not only avoid the contamination of the yu type distortions from the intermediate-type distortions 
but also measure the intermediate-type distortions themselves. 

The issue of degeneracy and the sensitivity required to detect intermediate type distortions can 
be made more precise by asking a slightly different question: How closely can a total spectrum 
containing all /i, y and intermediate-type distortions can be fitted by just fi and 3'-type distortions. We 
showinFig. 15 the total spectrum from Silk damping with «ss - 1 , and amplitude lOA^ = 1.61x10"^. 
The total spectrum has a Bose-Einstein part with fi = 10~^ and j-type part with y = 3 x 10"^. In 
reality, the y-type part of the spectrum would be much higher because of the contributions from 
reionization and later times, but this does not change our arguments or conclusions. Also shown is 
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Figure 15. Total spectrum containing an intermediate-type distortion from Silk damping with «ss = 1, and 
amplitude lOA^ = 1.61 x 10"*^, and jj and y-type distortions with ju - 10"^ and y = 3 x 10 All three 
components have approximately the similar amount of energy. Also shown is the best fit y^+ji spectrum to the 
data points at PIXIE frequencies. The best fit spectrum has = 1.8 x 10"^ and y = 4.26 x 10"**. The residual 
(data-fit) is shown in the bottom panel. 



the least squares fit using tiie data points at PIXIE frequencies (also shown) to the spectrum with only 
y and //-type distortions. The best fit spectrum has y = 4.26 x 10~^ and /i = 1.8 x 10"^. The bottom 
panel shows the difference between the data points and y + fi&t It is clear from these plots that the 
full spectrum cannot be fit exactly with only y and //-type components and the residuals (data-fit) 
contain give information about the type of intermediate-distortion present. In particular the residuals, 
which are ~ 20% of the intermediate-type distortion, are not affected by the presence of additional 
y and ju components in the spectrum. Thus the detection of the intermediate-type spectrum will be 
challenging but possible. In addition, since we do not know the temperature of the blackbody part of 
the spectrum at the required precision a priori, we should fit for the temperature of blackbody along 
with the distortions [1]. 



9 Observational issues 

Our definition of spectral distortions and xq is convenient to compare and understand the process of 
comptonization and theoretical spectra of different types. Observationally, it is not possible to calcu- 
late the reference blackbody temperature required for these definitions to desired accuracy. Achieving 
an accuracy of ~ 10~^ in the number density of photons, and hence the reference temperature, re- 
quires integrating the spectrum between 10^"^ < x < 25! This means that xq, Xminj -"^max defined 
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by us is not a direct observable for small distortions, xq, Xi-^m, ^max can however be inferred once the 
full spectrum is measured with sufficient precision and the intermediate-type spectral distortion part 
separated. 

The shape of the intermediate-type spectral distortions contains much more, and complemen- 
tary, information compared to the fi /v-type distortion. High precision experiments in the future may 
thus be able to not only put stringent constraints on the energy release in the early Universe us- 
ing /J. distortions, but also distinguish between different mechanisms of energy injection using the 
intermediate-type spectral distortions. Detecting the yu, y and intermediate-type distortions at the 
level of 10"^ - 10"'^ would require understanding and subtracting foregrounds at high precision. Sim- 
ulations for the proposed experiment PIXIE [5] show that such a subtraction may be possible at nK 
level. However, there aie still some uncertainties in our understanding of the foregrounds as demon- 
strated by the unexplained excess flux at low frequencies in the ARCADE experiment measurements 
[79] and more work needs to be done to demonstrate the feasibility of the measurements proposed in 
the present paper. 

10 Conclusions 

The CMB blackbody spectrum in the AC DM cosmology is established at z > 2 X 10^, when the 
photon number changing processes of bremsstrahlung and double Compton scattering are effective. 
Compton scattering establishes a Bose-Einstein spectrum corresponding to the energy and number 
density of photons available while bremsstrahlung and double Compton scattering help drive the 
chemical potential to zero. If there is energy/photon production due to non-standard physics at z < 2 x 
10^, the blackbody spectrum cannot be restored and an imprint is left on the CMB spectrum. Initially 
a j-type distortion is established which evolves (comptonizes) towai^ds a Bose-Einstein spectrum. 
How close the initial y-type distortion can get to a Bose-Einstein spectrum depends on the redshift 
of energy release. At z > 2 x 10^, full comptonization is possible but at lower redshifts Compton 
scattering can no longer establish a Bose-Einstein spectrum. At redshifts 1.5 x 10^ < z < 2 x 10^, 
there is partial comptonization and we get a spectrum which is in-between a j-type spectrum and a 
Bose-Einstein spectrum. The detailed shape of this intermediate-type spectrum depends on how the 
rate of energy injection varies with the redshift. The detection of a deviation of the CMB spectrum 
from a blackbody thus contains information about the amount as well as the redshift/mechanism of 
the energy release. 

We have numerically calculated the detailed evolution of the initial y-type distortion by solving 
the Kompaneets equation taking into account the con^ect evolution of the electron temperature, which 
decreases from the initial y-type distortion value towards a Bose-Einstein value. An analytic solution 
vaUd in the weak comptonization Umit, yy < 0.1, is given in Eq. (5.3). We have demonstrated 
that the intermediate-type spectral distortions, resulting from sound wave dissipation in the early 
Universe can, in principle, constrain the shape of the small scale power spectrum. It is also possible 
to distinguish between different energy injection mechanisms, which have different dependence on 
redshift, for example, WIMP annihilation and Silk damping with different power law dependence 
and particle decay with exponential dependence. The y, jj. and intermediate-type distortions have 
important differences in their shapes and it is in principle possible to distinguish between them. In 
particular-, a mixture of y and //-type distortions fails to mimic the intermediate-type distortions by 
20%, and vice versa. Proposed experiment PIXIE [5] will detect y and //-type distortions at the level 
of y = I0~^,iu = 5 X 10 *^ improving the current limits by 3 orders of magnitude. The additional 
y-distortions from low redshifts should not affect an experiment's ability to detect fj. and intermediate- 
type distortions. Measurement of intermediate-type distortions will thus be challenging but possible. 
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We make publicly available^ high precision intermediate-type distortion templates and a Math- 
ematica code which superposes these templates, according to user-defined redshift-dependent energy 
injection rate, to calculate the fi, y and intermediate-type distortions. 
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A Fitting formulae for xq, Xy^m, -^max as a function of yy 

xo(yy), .XminCVy), JCmaxCVy) fof the intermediate type spectra in Fig. 7 are well fitted (for small distor- 
tions) with an accuracy better than 1% for yy < 10 by the following formulae: 

XO/min/max = ^0 + aijy + ^2^7^ + «'33'r^ + '^^Jy'^ + ^1 sintT^bjyy) 

+ ci tanh(c23'y), (A.l) 

where we have for xq, 

ao = 3.83, ai = 0.363,a2 = -5.68 x 10~^^^3 = 5.15 x 10"^^^4 - -1-85 x 10""^, 

bi = 0.496, b2 = -30.8, ci - 0.294, C2 - 20.6, (A.2) 
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for JCmin 



c?o = 2.265, c?i = 0.332, £22 ^ -5.07 x 10"^ 03 = 4.4 x 10~\a4 = -1.5 x 10""^, 
bi = 0.439,^2 - -21, ci = 0.24, C2 - 14.1, 



(A.3) 



and for Xr, 



ao = 6.51, fli = 0.543,^2 = -0.11, as = 1.1 x 10"^£^4 ^ -4 x lO"'^, 
= 0.5, b2 = -89, ci - 0.56, C2 - 51. 



(A.4) 



B Analytic approximate solutions of Kompaneets equation 



We expand the photon occupation number (n(x, yy)) around the initial black body spectrum at tem- 
perature T, n{x,0) = Hpiix) = \/{e'"'""^'^ - 1) = l/{e^ - 1). 



n{x,yy) = n{x, 0) + yy—{x, 0) + ^^^C-^' 0) + 0(yy ), 



dn 



(B.l) 



where x = hv/ksT and T is the initial blackbody temperature. Further we can evaluate the corrections 
using Kompaneets equation. 
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(B.2) 



where ny{x) is just the linear j-type solution found in [3]. The next term in the Taylor series is 



dyy^ ' dyy yp dx 

-A - — 

' 7?- dx 



(^r, + l) + n + 

^^ dn{x,0) 1 Mr,^^^ 
dx At^ dyy 



yy=o 



+ At ——X 



1 d 4dMx,0) 



2 dx dx^ 



X ^ dx\ dx I \ dx I 
= ATjyix) + Aj^hix) + ^(0)nyix). 



(B.3) 



The functions fy{x) and f2(x) are given explicitly in the Appendix C. Substituting the first and second 
order terms from Eqs. (B.3) and (3.4) in Eq. (B.l), we have the solution for partial comptonization 
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of an initial blackbody spectrum interacting with electrons at temperature T^iyy), valid for small 
distortions and correct to second order in yy, 



n{x,yy) =?ipi +yy^Tny 

2 



+ 



yy 



dyy 



+ 0(yy\ 



(B.4) 



where Aj^ and its derivative are evaluated at yy = 0. 

The last term, ^p-, depends on the physics responsible for the electron temperature Tg. It is 
also clear that higher order terms will have higher order derivatives of Aj^ and the above solution is 
only valid when these higher order derivatives are negligible compared to the first derivative. We will 
ignore this term in the rest of this section, as it is not relevant to the present discussion. In the early 
Universe, however, electron temperature is not constant and changes as comptonization progresses 
and we will discuss the evolution of the electron temperature in detail in the next sections. If the 
electrons are in equilibrium with radiation, is easily calculated by taking derivative of Eq. (5.1). 

We can now use the fact that the blackbody spectrum is a steady state solution of the Kompa- 
neets equation. 

dn{x, 0) 



n(x, 0) + n(x, Oy = 



dx 
dn(x, yy) 



^ " . n^ yy 



d^n 



■ix,0) 



dx ^ '' dyydx ^ ' ' 2 dyy^dx 
+ 0(yy') (B.5) 

We have used Eq. B.l in the last step. Using again the expansion Eq. B.l and Eq. B.5 in the 
Kompaneets equation to replace n + rP- term and , we get, on ignoring terms of order yy^ and higher, 

dn{x,yy) 1 d 



dyy 



x^ dx 
+yy 



dn{x,yy) /Te _ \ 
dx \T ) 



d^n dn dn 

— — (X, 0) + —{X, 0) + 2n{x, 0)—{x, 0) 
dyydx dyy dyy 

d^n d^n d^n 

(x, 0) + — ^(x, 0) + 2n{x, 0)— ^(x, 0) 



2 \ dyy^dx 

I dn 
+2\—(x,0) 



dyy^ 



'dyy^ 



(B.6) 



We are interested in the behavior of corrections as yy and Aj^ increase. Thus collecting all terms of 
same order in Vy and Aj^ (and defining representing x-dependence with functions F„{x) for brevity)^, 
we have 



dn{x,yy) _ \ d 
x^ dx 



Ar.dy 



y 



^dn{x,yy) 
dx 



+ 0(yy') 



x»i 1 d 
x^ dx 



^ dn{x,yy) 
dx 



+ yyFi{x) + yy'^F2{x) + yy^AjF^ix) 



{lyyX^ - "iyy^x^ + lyy^ Aj^x'^n.g\{x) 



+ (5(3'r ) 



(B.7) 



(B.8) 



*We only give the x » \im\1 below but the functions F„(x) just involve the Planck function ;7pi and its derivatives are 
easily calculated explicitly if desired. 
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The last expression illustrates clearly the regime of validity of the y-type solution. The y-type solution 
is valid in the limit <sc 1 and y~f-lS.j^ «c 1. Also the solution fails at ;c » 1, when the recoil effect 
gives deviations of order unity with respect to the blackbody. The recoil effect would lead to a 
downward shift in the frequency of high energy photons given \)y \ jx' - \ jx = yy [80, 81]. Thus the 
j-type solution is only valid for x «; 1 jyy. This is also easily seen by comparing the first term on the 
right hand side (~ x^?ipi) with the first two terms in round brackets, Q.yyX - 3yy^x^)x^«pi in the limit 
X » 1. 

In particular for A^^ «: 1 , as is the case in the early Universe before recombination, the «: 1 
is the stronger condition. Thus the 3'-type solution is the correct solution for energy injection at 
z < 20000 with the connections due to higher order terms of the order ~ yy ~ 10"^ (Fig. 4). 

With y - t^j^Ayy and y 1, an approximate linear solution of Eq. (B.7) follows by 
evaluating the right hand side of Eq. (B.7) at yy = 0, with y = I /4{AE/Er) , AE is the energy injected 
into the plasma and E,- is the initial radiation density. 

nix,y) - ?ipi(x) = yuyix) 

xe^ r e-^ + 1 1 
= y T X --4 (B.9) 

An important difference between Eq. (B.9) and linear + second order AAjJAyy term in Eq. 
(B.4), although they look identical, is that in Eq. (B.9) we have defined y as an integral over A^^ and 
there is no restriction on the higher order derivatives of Aj-^ , except that the integral y «c 1 and the 
assumptions under which Kompaneets equation is derived are valid. We have in effect summed over 
all the terms coming from the Taylor series expansion of A^^ in Eq. (B.4). The solution in Eq. (B.9), 
written in terms of y, is thus, more generally applicable. This solution is valid for An/upi «: 1, which 
implies that for large x we have the condition x^y <«c 1. This condition is clearly satisfied for the 
CMB, with the current limit of y < 10"^, in the Wien tail for x < 100. The well known solution Eq. 
(B.9) [3] depends only on y and not on yy. The next correction depends also on yy. This demonstrates 
that the broadening of the spectrum and redistribution of photons over the frequency is defined by 
but the energy exchange between the plasma and the radiation is defined by Tg - T. 



C Corrections to y-type distortion from weak comptonization and recursion relations 
for calculating the higher order terms 

The solution to the Kompaneets equation for small distortions, with the initial spectrum being a 
blackbody spectrum, for small values of yy parameter, yy «; 1, is given by the Taylor series Eq. 
(B.l), with the first two coefficients/derivatives given by Eqs. (3.4) and (B.3). The functions fyix) 
and f2(x) just involve derivatives of the Planck spectrum and are easily calculated. 

f^,{x) = \-x^ -Ix-S + e^' (x^ -Ix+S) 

ie' - 1)5 L ^ ^ 

[-2x^ - 9x^ + lx + 24) + e^"" [-2x^ + 9x^ +lx- 24)] (C. 1) 

fiix) = ^ ^ , \x^ + 12^2 + 34x + 16 + e^'^ix^ - I2x^ + 34x - 16] 

+£-2' (l Ix^ - 36x2 _ ^ 43 j ^ Jhx^ + 36x2 _ 34^ _ 43 jj ^^2) 

For an initial y-type spectrum, ?i(x, 0) = «pi(x) +y?ij,(x), Aj^iyy = 0) = 5.4y, we can similarly 
expand the solution n{x,yy) in Taylor series, keeping only terms linear in y since COBE/FIRAS [1] 
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100 1000 




Figure 16. Analytic solution given by Eq. (C.4) including terms up to is shown by dotted lines for 
yy = 0.0 1 , 0.05 , 0. 1 from top to bottom (in the Rayleigh-Jeans and Wien tails) respectively. Numerical solutions 
are as marked. Fractional difference in the effective temperature, Eq. (3.5), is plotted. The two solutions match 
well forjy «; 1 and x «; 1 /yy. At yy = 0.1 the analytic solutions has the correct approximate shape. Errors are 
plotted in Fig. 17 



already constrains the average cosmological distortion, y < 10"^, yy on the other hand covers a wide 
range and is > 1 at z > 1.45 x 10^). yy in the solutions refers to the total yy integrated from the time 
of energy injection to the time where we want to calculate the final distortion. 



dn 

' dyy 



n{x,yy) = Upiix) + yuyix) + yy—ix,0) + + + ^^Vr ^ 



(C.3) 



We can calculate the coefficients in Taylor series iteratively by using Kompaneets equation. The 
first derivative is thus given by the Kompaneets equation. Second derivative is obtained by dif- 
ferentiating Kompaneets equation with respect to Vy and using the solution of first derivative and 
so on. Derivatives of Aj^iyy) are also easily obtained using Eq. (5.1). First two derivatives are 
dArJdyyly^^o ~ -21.45^, and <fATjdyy\=o « 323.6);. Thus, 



n{x,yy) = npi+y 



y + yy (5 Any +fy) + Y' (-21-45«y + 5.4^. + ) 



+ZZ_ (323.6«, - 21.45/, + 5.44^) + ^3)^ 



+ 0(V). 



(C.4) 
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Figure 17. Error in analytic solution defined by ^^'^^^^Tm ^'^^^""""'"'■^ . The analytic solution has better than 
1% accuracy over most of the frequency range of interest at yy < 0.05. Accuracy deteriorates quickly at larger 
values of yy. 

In general derivative is given by (for m > 1), 



where gf^ - riy, gyP - fy and for m > the gy'"'' functions are given recursively by 



(m+l) _ 1 4 



— +g^"'\l+2«pi) . (C.6) 



dx 



We compare the the analytic solution including first three terms (up to order yy^) and the nu- 
merical solution in Fig. 16. The two solutions match very well for yy ^ I, x <^ l/yy and the error 
at yy = 0.05 is < 1% for x < 7. However, for larger values of yy, the solution quickly deteriorates, 
and at yy = 0.1 the error is of order 10% at ;c < 6. For yy < 0.01, the linear order term is enough 
to give ~ 1% accuracy. Fig. 18 compares the analytic solution including up to linear, quadratic and 
cubic terms with numerical solution for yy = 0.1. It can be seen that with the inclusion of successive 
terms, the analytic solution oscillates around the true solution and convergence is quite slow. Table 
1 gives the approximate maximum values of yy where the error at < 6 is below 1%,5%, 10% at 
linear, quadratic and cubic orders in yy. 
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Figure 18. Error in analytic solution defined by ,^T!^, °Tn"^^ for - 0.1 for analytic solutions at 

different orders. The analytic solution has better than 10% at x < 6. It can be seen that the convergence of 
the Taylor series is very slow at = 0.1 with the analytic solution oscillating around the true solution with 
inclusion of successive terms. 



Error 


1% 


5% 


10% 


linear 


0.01 


0.03 


0.04 


quadratic 


0.03 


0.06 


0.08 


cubic 


0.05 


0.08 


0.1 



Table 1. Upper bound on errors in the analytic solutions at different orders in yy. The upper bounds in different 
columns are reached at the values of shown. This table can be used as a rough guide when using analytic 
solutions. 
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